Circular Optical Nanoantennas - An Analytical Theory 



(N. 

o- 
<n: 

a- 
»— ». 



o- 

-<— > 

q. 

>^: 

Or 



>: 

av 

o: 



X: 




Robert Filter, Jing Qi (M$W), Carsten Rockstuhl, and Falk Lederer 

Institute of Condensed Matter Theory and Solid State Optics, Abbe Center of Photonics, 

Friedrich-Schiller-Universitat Jena, Max-Wien-Platz 1, D-07743 Jena, Germany 

An analytical approach is provided for describing the resonance properties of optical nanoantennas 
made of a stack of homogeneous discs, i.e. circular patch nanoantennas. It consists in analytically 
calculating the phase accumulation of surface plasmon polaritons across the resonator and an addi- 
tional contribution from the complex reflection coefficient at the antenna termination. This makes 
the theory self-contained with no need for fitting parameters. The very antenna resonances are then 
explained by a simple Fabry-Perot resonator model. Predictions are compared to rigorous simula- 
tions and show excellent agreement. Using this analytical model, circular antennas can be tuned by 
varying the composition of the stack. 
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I. INTRODUCTION 

Ever since Hertz described in 1887 the emission of elec- 
tromagnetic radiation from dipole antennas, a pertinent 
question has been to fabricate them such that they may 
interact with light at optical frequencies, a spectral do- 
main of paramount importance for many applications. 
To achieve this goal antennas have to be downscaled in 
their critical dimensions to a few hundreds of nanometers 
to match the wavelength of the visible. More than a cen- 
tury later, this spectral domain has been reached owing 
to recent advances in nanofabrication and characteriza- 
tion techniques b However, despite of this progress, the 
theoretical understanding of optical nanoantennas is lag- 
ging behind the available technology. The complexity of 
this issue is also one of the reasons why the topic attracts 
much research interest 2 . Whereas metallic antennas at 
radio frequencies, where the metal may be considered as 
perfect conductor, can be conveniently treated analyti- 
cally using Babinet's principle, this ceases to hold in the 
visible 3- -. There, the dielectric functions of metals have 
to be properly accounted for and may be phenomenolog- 
ically described by the free-electron (plasma) model. 

The arising plasma oscillations may couple to the elec- 
tromagnetic radiation at a metallic interface to form 
new confined quasi-particles 6,7 . These quasi-particles 
are termed surface plasmon polaritons (SPPs) and gov- 
ern the resonant behavior of optical nanoantennas 8 . 
Such optical nanoantennas enabled already ubiquitous 
applications^— and will certainly find many future ones. 

Eigenmodes of various nanoantennas have been mea- 
sured by different meansi^— . However, problematic in 
further advancing the field is the lack of analytical insight 
into the scaling behavior of optical nanoantennas to care- 
fully design them for a desired application 2 ^—. Most de- 
sign related tasks rely entirely on numerical tools; hence 
providing only little insight into the underlying physics. 
This is in stark contrast to the field of antennas at ra- 
dio frequencies and it would be highly desirable to bridge 
this gap. 

First efforts in this direction has been done character- 



FIG. f. The excitation of a Bessel-type SPP at a specific 
circular nanoantenna, a metallic disc. The reflection at the 
termination leads to an apparent length change as described 
by the given theory. 



izing the resonance behavior of the potentially simplest 
optical nanoantenna, i.e. the nanowire antenna. It was 
shown that their resonances can be explained by a Fabry- 
Perot model. Envisaging the nanowire as a plasmonic 
cavity, a resonance occurs whenever the phase accumu- 
lated upon a single round trip amounts to a multiple of 
2-7T. This phase is determined by the propagation con- 
stant of the plasmonic mode supported by the nanowire 
and an additional phase upon reflection at the wire ter- 
mination. By accounting for that additional phase, the 
length of the nanoantenna as perceived optically differs 
from its geometrical length. This phenomenon is some- 
times termed the apparent length change of the nanoan- 
tenna which turned out to be essential since it strongly 
affects the resonance positions 3 ^—. In the past this ad- 
ditional phase contribution has been considered either ad 
hoc without further justification or as a free parameter in 
an adapted model. By fitting many resonance positions, 
taken from many samples with different nanowire lengths 
and/or incidence angles, to predictions from the model, 
the phase change upon reflection could be phenomeno- 
logically determined. Such tedious procedures were nec- 
essary since the analytical calculation of the complex re- 
flection coefficient turned out to be fairly involved. To 
date such results are only available for a few simplified 
geometrie s 36 ' 38 . Understanding and predicting the com- 





FIG. 3. A Hankel-type SPP propagates outwards a circu- 
lar nanoantenna and gets reflected at its termination into an 
inwards propagating Hankel-type SPP. 



FIG. 2. Schematic of a circular nanoresonator. Its permit- 
tivity e (z) is axially symmetric and piecewise homogeneous. 
The whole structure is radially terminated at p — R and em- 
bedded in a medium of permittivity td- 



plex reflection coefficient of surface plasmons at abrupt 
interfaces is one of the key issues to be solved for a future 
advance in nanophotonics. 

Our contribution towards this goal is twofold in con- 
sidering here circular patch nanoantennas, i.e. nanoan- 
tcnnas made from an arbitrary stack of homogeneous 
discs. First, the complex reflection coefficient of Hankel- 
type surface plasmon polaritons at the circumference 
of the circular patch nanoantenna is analytically calcu- 
lated. This consideration is extremely versatile since 
many other, but simpler geometries, can be analyzed 
while probing for limiting cases of this theory. For ex- 
ample, increasing the radius of a single disc towards in- 
finity and increasing its thickness yields the reflection 
coefficient of an SPP propagating along a single metallic 
interface at a planar termination. 

Moreover, this theory can be used to predict and engi- 
neer the resonances of circular patch nanoantennas made 
of an, in principle, arbitrary sequence of layers. This is 
achieved by a Bessel-type standing wave resonator model 
using the calculated phase of reflection (see Fig. Q] for a 
graphical illustration). Analytical results will be com- 
pared to both known results and rigorous simulations to 
prove the applicability of this theory. 

With this work a simple analytical tool is provided 
that allows for a deeper insight into the scaling behavior 
of complex nanoantennas and, moreover, defines a path 
how to treat such nanooptical systems analytically where 
only numerical approaches have been ruling in the past. 



II. REFLECTION OF HANKEL-TYPE SPPS 

The class of nanoantennas discussed in this work is il- 
lustrated in Fig. [2] Such axially symmetric structures 
support radially propagating SPPs which are described 
by Hankel functions and consequentially termed Hankel- 
type SPPs^. In this section, the reflection coefficient 



of such SPPs at the circumference of axially symmetric 
nanoantennas will be introduced. Because of the singu- 
larity of the Hankel functions at the origin, this ansatz 
cannot be used to describe the actual modes of circular 
nanoantennas which will be discussed in the following 
section. 

Assuming a linear response of all involved materials 
it is sufficient to consider time-harmonic solutions with 
a exp (— iuit) time-dependence. Furthermore, the plas- 
monic modes obey an angular variation of exp (irrup) 
where m denotes the angular mode index. In the fol- 
lowing, the time dependency will be dropped and the 
angular dependency will be denoted by the subscript m. 

Assuming a transverse-magnetic (TM, H z = 0) field, 
it is sufficient to describe the fields in terms of E z alone. 
The plasmonic modes obey the usual boundary condi- 
tions at discontinuities of e (z) and the propagation con- 
stant in radial direction is subject to a dispersion relation 
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For the given time-dependency, the Hankel functions 

1/2 

Hm of the first and second kind correspond to outwards 
and inwards propagating cylindrical waves. For Hankel- 
type SPPs, the radially propagating parts of E z are then 

1/2 

described by a superposition of the H„[ . 

Outwards propagating modes are reflected into inwards 
propagating ones at the circumference of the nanoanan- 
tenna. The approach for calculating the corresponding 
reflection coefficients r m along with the results will be 
briefly presented in the following. The detailed deriva- 
tion can be found in App.[SJ A schematic of the reflection 
is given in Fig. [3J 

For the described reflection, the radial part of the sur- 
face plasmon field can be described by 
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for p < R implied by the superscript — ; + will be used 
later to denote p > R, the region outside the nanoan- 
tenna. The mode profile a (z) determines the propagation 
constant kspp. Furthermore, It is sufficient to calculate 
r m for a fixed angular mode m since there will be no re- 
flection into modes with a different rotational symmetry. 



It must be noted that this description implicitly ne- 
glects the reflection into other eigenfunctions of the res- 
onator and that for strong absorption this Hankel-type 
form ceases to be valid; being essentially replaced by Nor- 
ton waves as recently describe d 40 ' 41 and further discussed 
in section HVBI 

Outside the resonator the field couples, as in the case 
of reflection, only to modes of the same angular symme- 
try class. A suitable representation is a continuum of 
outward propagating Hankel waves 



E'l 



(r) = 



[k z )H x m [Je d kl-klp )e?"*dk 



with the corresponding expressions for H™ ,+ . The am- 
plitudes c m (k z ) are given by obeying the boundary con- 
ditions and any scattering of objects outside the nanoan- 
tenna e.g. by Bragg mirrors 4 ^, is neglected. 

Following Gordon^, the reflection coefficients r m are 
determined using the continuity of H v and J E z ■ H v dz 
at p = R. Thus, the complex reflection coefficients read 
as: 

_ 2ne d ksppaH^jksppR) - DH^(k SPP R) I m 
Tm -2ire d k S ppaH? n {k S ppR)+DHl{k S ppR)I m U 

with the abbreviations 
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To compute the phase change upon reflection <$> r m , the 
origin of the coordinate system has to be specified. As 
a result the phase directly coiuputed from r m is not 
the phase upon reflection, but an overall accumulated 
phase </>^ c with respect to the origin of the coordinate 
system. This overall phase can be decomposed into a 
reflection (jf m and a propagation contribution <^, thus 
r m = \r m \ exp (i^T) with ^T = C + C • For Hankel- 
type fields, one finds 



= Arg [H} n (fc S pp • 0+)* • H x m (fcgpp • R) 
Arg [Hi (fc S pp • R)* ■ H 2 m (fc SPP • 0+)] 



where a further decomposition of the phase accumula- 
tion into in- and outwards propagating waves was used. 
Equation [5] is the main analytical result of this work and 
constitutes the desired analytical expression for the com- 
plex reflection coefficient. 



III. A SIMPLE RESONATOR MODEL 

To calculate the complex reflection coefficient it was 
necessary to use Hankel-type SPPs. Otherwise, in- and 
outwards propagating waves could not be defined. Un- 
fortunately, these functions are not finite at the origin 
and thus not meaningful as standing-wave solutions for 
E z inside the resonator. Thus, to describe the actual 
field one has to use a plasmonic mode that takes finite 
values at the origin. These modes are Bessel-type SPPs 
described in analogy to Eq. Q] via 
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(p, z) = J m (fcsppp) • a (z) 
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The question is how one can combine these two descrip- 
tions in the sense of interpreting propagating SPPs as lo- 
calized ones 4 ^,. As noted earlier, the phase change upon 
reflection of one-dimensional nanoantennas causes a de- 
viation of the resonance condition. The apparent length 
change was explained by the reflection phase and the 
same approach will be applied here. For a reflection with- 
out phase accumulation, the resonant radii R n ^ m of order 
n would be given by the roots of E™'~ leading to the 
roots of Bessel functions as for circular microstrip patch 
antennas 44 . Hence, introducing a nonzero phase upon 
reflection, one may conjecture the real- valued resonance 
condition 



Z • kgppR nrn + (f> rn — 2 ■ x n \J m ) 
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where x n denotes the n-th root of J m , fcspp = fcg PP + 
ifcgpp. The given formulation is a natural extension of 
the Fabry-Perot resonance condition for one-dimensional 
nanowire antennas 35 . There, the SPP's were character- 
ized by guided plasmonic waves confined by the nanowire 
and whose propagation along the nanowire is character- 
ized by a propagation constant fcspp- Resonances were 
found by solving 2fcg PP L„ + 2(/> r = 2im. It must be noted 
that for reflection amplitudes that strongly deviate from 
unity, a feasible description of the excitation in terms 
of localized plasmonic modes might not be possible any- 
more. Thus the employed resonator model loses its pre- 
dictive strength. 



IV. VERIFICATION OF THE THEORY 

The latter two equations for the actual form of the 
fields and for the resonances of circular patch nanoan- 
tennas are just educated guesses at this point that have 
to be verified. 

Although different numerical methods could have been 
used to achieve this 4 ^— , finite-difference time-domain 
(FDTD) simulation s 48 ' 49 were carried out for the follow- 
ing subsections. The obtained results were compared to 
theoretical predictions. 

In the first subsection, the assumed Bessel form and 
scaling of E z is studied and a comparison to the pre- 
dictions for resonant radii is performed. Metallic discs of 



different thicknesses are used as circular nanoantennas at 
a fixed frequency. A very good agreement is shown. The 
resonator model can further be used to calculate the res- 
onance frequencies of a circular nanoantenna with fixed 
geometrical parameters. An agreement to simulations is 
also found in this case and outlined in the second subsec- 
tion. Furthermore, it is shown how derived quantities can 
be used to calculate the quality factor (Q-factor) of the 
nanoantenna. Limitations of the theory are discussed. 
The section will finish with a study of a stacked structure 
consisting of two silver discs with a cavity in-between. It 
will be shown that the theory is also able to accurately 
describe this situation. This subsection will prove that 
the theory is versatile and that it can be applied also to 
circular nanoantennas with a larger complexity. 

The convergence to known results in limiting cases is 
outlined in App. [Bj 



A. Eigenfunctions and Scaling 

Equations [3] and 0] imply two different things. First of 
all, the radial dependence E z is assumed to be described 
by a Bessel function with a certain scaling caused by the 
dispersion relation. Second, resonant radii are predicted 
for which the plasmonic field enhancement is strongest. 
Both assumptions will be tested in this subsection at a 
fixed frequency where the permittivity of the metal is 
just a constant and the observed effects are independent 
of any nonlinear spectral scaling of fcgpp(w). 

Discs of several thicknesses di and radii between 
100 nm and 1250 nm with a relative permittivity e m ~ 
— 8.8 4- 0.03i were illuminated by a plane wave propagat- 
ing normally to the disc surface and having a frequency 
of v = 625 THz. The material parameters correspond 
to those of silver at this frequency 50 . The imaginary 
part has been deliberately lowered to simplify the iden- 
tification of the resonances. However, for all subsequent 
subsections the values were directly taken from the liter- 
ature. Because of the known wave vector mismatch, the 
plasmonic modes on planar discs do not couple to those of 
free space. Hence, the excitation of the SPPs generally 
takes place at the disc's edges. Furthermore, a certain 
mode can only be excited if the outer field has approx- 
imately the same spatial symmetry where the coupling 
strength is largest. Although the discs are sufficiently 
thin such that in general even and odd propagating sur- 
face plasmons are supported, the latter ones could not 
be excited since the disc thicknesses were much smaller 
than the wavelength. The required antisymmetric field 
distribution along the discs thickness does not match to 
the incident plane wave which shows a phase variation on 
length scales in the order of the wavelength. Hence, only 
even SPPs could be excited with a reasonable efficiency. 

The resonance radii R n = i?„,i were determined by 
identifying maxima of the electric field strengths below 
and above the structure while changing the radius of discs 
at a constant thickness. The verification of the employed 
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FIG. 4. Numerical verification of the circular Fabry-Perot 
resonator model, a) The electric field component E z on a 
metallic disc at resonance shows a qualitative agreement with 
the assumed Bessel-type form, b) Resonant radii R„ for disc 
thicknesses of 6 nm, 10 nm, 20 nm, 40 nm, 80 nm, and 160 nm 
are linearly related to the roots of Ji. c) The normalized 
field E z on a metallic disc for the first (blue line) and fourth 
(magenta line) resonant radius. Inside the discs p < R the 
agreement with the expected form oc J\ (ksppp) (green dashed 
line) is perfect except a minor deviation at the termination, 
fcspp was calculated from the known dispersion relation for a 
metallic slab&i. 



resonator model requires that a) the excited fields exhibit 
the assumed form and b) the measured R n are linearly 
related to the roots of J\ . Looking at Fig. QJ it can be rec- 
ognized that both conditions are met. There, the fields as 
predicted by Eq. [3] are compared to full wave simulations 
and an excellent agreement is observed. The fields at 
nanoantenna resonances follow exactly such Bessel-type 
SPPs. Moreover, the resonance radii are linearly related 
to the roots of J\ as assumed in Eq. |U Consequently, 
Eqs. [3] and U seem to provide an appropriate description 
of the actual physical situation. 

The most interesting question might be if the resonant 
radii R n can be explained by an analytically calculated 
phase of the reflection coefficient for the excited even 
modes by using Eq. 0] and thus the phases of the com- 
plex reflection coefficient. The calculated phases <p\ can 
be seen in Fig. [5] Phases corresponding to different thick- 
nesses are significantly different. For the given results, it 
seems to be a monotonic behavior between thickness d 
and phase upon reflection <p\ - the thicker the disc the 
larger the reflection phase, </>[ (di) < <p\ {d^) < • • ■ for 
di < dj+i. Furthermore, all c\>\ turn out to be negative, 
the only exception arises by the two thickest discs for 
small radii. Noteworthy, the relative change of cf>[ from 
the first to the fifth resonance radius amounts approxi- 
mately 10% to 30%. It is increasing with thickness. 

The comparison of numerically and analytically calcu- 
lated resonant radii is displayed in Fig. [51 The predictions 
for even modes coincide very nicely with numerical sim- 
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FIG. 5. Phases upon reflection vs. the disc radius calculated 
from the theoretical model for metallic discs of the same thick- 
nesses as in Fig. [4] The </>i turn out almost always negative. 
With respect to the existing terminology 35 , this corresponds 
to an apparent length decrease. The effect increases for thin- 
ner discs. 




20 40 80 160 



(Kinii) 



FIG. 6. Analytically calculated phases of the reflection co- 
efficient are used to predict the resonance radii R n employ- 
ing Eq. |4] (full lines) for discs of different thicknesses d. The 
five lowest order resonances are shown. Results from rigorous 
simulations are shown by dots. Large squares correspond to 
R n resulting from calculations assuming a semi-infinite disc; 
hence the thickness d is infinite. 



ulations for thicknesses up to 80 nm. This suggests that 
for a given thickness the theory predicts a series of reso- 
nance radii which are in full agreement with those found 
using numerical simulations. For an increasing thickness 
one can observe that the resonance radii are independent 
of the thickness. This can be explained while considering 
the characteristic penetration depths of the fields into the 
metallic discs which are in the order of 20 nm. Evidently, 
for thicker discs the fields confined at the top and the 
bottom interface are decoupled. Thus the reflection co- 
efficients coincide then with those of semi-infinite discs. 
This explains why for thicknesses exceeding 40 nm the 
resonant radii agree very well with the predictions of this 
approximation (see App. [Dlfor an explicit calculation of 
r m in this case). 



B. Spectral Dependence of Resonance Frequencies 



In the last subsection the comparison of theoretical 
predictions to results from simulations was discussed at 
a fixed frequency to avoid an inevitable nonlinear scal- 
ing due to the frequency dependence of fcspp (w). So, 
it is natural to ask if also the spectral resonant behav- 
ior of the circular nanoantennas is in agreement with the 
given model considering a fixed geometry. Furthermore, 
using the permittivity of a real metal, one can question 
what happens if reasonable losses at high frequencies are 
properly considered. Then, the theory should cease to 
be valid since mode profiles will differ from the assumed 
Bessel-type SPP. 

To answer these two questions, rigorous simulations 
were performed for silver discs where the material data 
was taken form literature^. The discs were assumed 
to have a radius of R = 900 nm and a thickness of 
d = 20 nm. They were illuminated by plane waves at 
different frequencies. Figure [7] a) shows the field max- 
ima above the structure as calculated from FDTD 48 i 49 
simulations as a function of the frequency. In addition 
the predicted resonance frequencies from the analytical 
theory are shown. An excellent agreement for nearly all 
orders can be noticed; although the model is slightly less 
accurate for higher order resonances. With respect to 
these results one finds first deviations for the fourth or- 
der mode at v k, 660 THz. There, the permittivity of 
silver is e m (660 THz) w —7.3 + 0.2 i leading to a wave 
vector of fc SPP « (1.6 • 10 7 + 1.2 • 10 5 i) m" 1 for even 
modes. The propagation length dspp = 1/fcgpp is in this 
case approximately 8.6 /jm which is roughly five times 
the diameter of the disc. Thus, this damping is not neg- 
ligible but with respect to the dimensions of the disc, the 
resonator model can be applied. 

The situation changes at the fifth order resonance at 
v w 740 THz. There, e m (740 THz) « -4.6 + 2.1 i and 
fc S pp « (2.0 • 10 7 + 3.5 • 10 5 i) m _1 leading to a prop- 
agation length of dspp ~ 2.8 /im which is in the order 
of the disc diameter. In this highly damped case, the 
resonance differs from the assumed Bessel-type SPP and 
the theory is no longer valid 4 ^ 4 ^. This can be also seen 
comparing the field amplitudes across the surface of the 
disc as shown in Figs. [7] b) and [7J c) for the fourth and 
fifth order mode. 

Looking at the resonance spectrum in Fig. 0a), one can 
pose the question if the theory can predict the width Av 
of the resonance peaks. Interpreting the Bessel-type SPP 
as damped oscillation across the nanoantenna which has 
been done for other kinds of resonant nanostructures 51 , 
this question directly leads to the quality factor Q5S. To 
calculate Q, the energy loss of the plasmonic mode per 
cycle has to be computed. Two direct contributions are 
given by the theory, propagation losses due to the imag- 
inary part of fcspp and radiation losses upon reflection. 
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FIG. 7. Comparison of theoretical predictions to simulation 
results for the resonance frequencies of a silver disc with ra- 
dius R = 900 nm and thickness d — 20 nm. a) The maximum 
value of \E Z \ as calculated from rigorous simulations as a 
function of the illumination frequency v (blue line). Reso- 
nances (dashed green lines) and line widths (magenta solid 
line) as predicted by the theory are additionally shown. For 
the first four resonances, one can observe a v^ scaling (yellow 
dashed lined) . This behavior very likely corresponds to a vari- 
ation of the incoupling efficiency which is beyond the scope of 
this work. Nevertheless it is worth mentioning that the scal- 
ing does not hold anymore for the fifth resonance and may 
be explained by the high damping. It furthermore indicates a 
need to adjust the theoretical description of the system in the 
respective regime. This is underlined by the actual plasmonic 
fields: b) The field distribution at tp — for the fourth or- 
der resonance at v ~ 660 THz (blue line) obeys the assumed 
Bessel-type form (green dashed line). In c) the situation is 
shown for the fifth order resonance at v ~ 740 THz. One can 
observe clear differences to the assumed form which can be 
attributed to a larger damping of the Bessel-type SPP. 



One finds 
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(2k'^ P Rf + 
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(5) 



Using Eq. [5J the quality factor for the Bessel-type SPP 
was calculated. The results are displayed in Fig. [8] Up to 
the fourth resonance, the quality factor is dominated by 
the radiative contribution Q T and only at higher frequen- 
cies damping of the Bessel-type SPP prevails and limits 
the quality factor. The predictions for the quality factor 
were used to calculate the peak widths Aza They agree 
very well with those seen in Fig. [7] a) for the numerically 
simulated response. 
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FIG. 8. The quality factor Q as calculated from Eq. [5] for a 
silver disc with radius R — 900 nm and thickness d = 40 nm. 
Up to 600 THz, losses due to propagation are much lower 
than radiative ones such that Q is entirely dominated by the 
contribution of Q T (magenta line). The reflection coefficient 
r is increasing with frequency in this spectral domain. Thus, 
Q T = |1 — r\~ shows the same behavior leading to narrowed 
peak widths as observed in Fig. [7] a). 



C. A Simple Stack System 

Until now for the sake of simplicity the theory was 
tested for an isolated metallic disc. However, the theory 
is also applicable for stacked nanoantennas as shown in 
Fig.[5J Such stacks introduce a much larger degree of free- 
dom for tailoring the properties of Bessel-type SPPs con- 
cerning both their propagation constant and their com- 
plex reflection coefficient at the disc termination. In this 
subsection, it will be shown that the theory is indeed ap- 
plicable to these kinds of systems by verifying it for a 
selected case. 

Namely, the nanoantenna under consideration shall 
consist of two silver discs, each of thickness d. They 
are separated by a dielectric spacer of the same thick- 
ness. For reasons of simplicity air was assumed as the 
spacer and the surrounding medium. The configuration 
is illustrated in Fig. [5] a). 

Generally the structure supports several Bessel-type 
plasmonic modes. However, only two modes can be ex- 
cited efficiently. They can be understood as the sym- 
metric and antisymmetric coupling of the metallic discs. 
Thus, the modes are termed symmetric and antisym- 
metric in correspondence to the symmetry class of E p 
with respect to the z-coordinate. If E p has a certain 
symmetry E z belongs to the opposite symmetry class - 
E z is symmetric if E p is antisymmetric and vice versa. 
Hence, E z vanishes for the symmetric mode in the center 
of the air cavity but is strongest there for the antisym- 
metric one. Furthermore, the symmetric mode is con- 
fined mostly at the lower and upper terminations of the 
structure. This spatial separation of the modes make it 
possible to identify resonances of the nanoantenna corre- 
sponding to the different modes. More on the nomencla- 
ture, dispersion relation and profile of the modes can be 
found in App. IC 31 
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FIG. 9. a) Schematic of the stacked nanoantenna - two sil- 
ver discs with an air cavity in-between. Every layer has a 
thickness d. b) For a non-resonant radius of R — 810 nm and 
thickness d = 20 nm, E z is plotted for the symmetric and anti- 
symmetric modes directly from simulations. The field profile 
for the symmetric mode on the left (orange frame) was com- 
puted 20 nm above the structure, whereas the antisymmetric 
one (blue frame) directly in the middle. Also, a comparison of 
the theoretical predictions (full lines) against simulation re- 
sults (dots) is shown. Again, a very good agreement is found. 
It has to be mentioned that for a thickness of d — 40 nm the 
resonances of the symmetric mode could not be identified. 
There, the antisymmetric mode turns out to be dominating. 



The simulations were analogue to the ones outlined in 
subsection IIV Al The nanoantennas were illuminated by 
a plane wave propagating in z-direction at a constant fre- 
quency v — 625 THz. There, the permittivity of silver- 
is e m ps —8.8 + 0.3L The radii R of the antennas varied 
from 40 nm up to 1000 nm and d was chosen to range 
from 6 nm to 40 nm. Absolute values of E z were moni- 
tored in the middle of the structure and 20 nm above and 
below to identify resonances of the odd and even modes, 
respectively. The results are outlined in Fig. |H] b) . 

Hence, it may be concluded that the reflection coeffi- 
cient as calculated from Eq. [5] can be used to evaluate 
the resonant behavior for the given class of nanoanten- 
nas if the dominant modes inside the resonator are the 



assumed Bessel-type SPPs. 



V. CONCLUSION 

A theory to describe radially propagating Hankel-type 
SPPs in piecewise homogeneous circular nanoresonators 
was introduced. The complex reflection coefficients at 
the termination were analytically calculated, thus yield- 
ing the phase change of the SPPs upon reflection. In 
combination with a Fabry-Perot model that predicts the 
resonances of circular patch nanoantennas, the theory 
was proven to be valid while comparing results to full- 
wave simulations. 

This theory constitutes a unique approach towards the 
analytical discussion of the resonant properties of circu- 
lar optical nanoantennas without requiring any fitting pa- 
rameters. This allows for a deeper insight into the scaling 
behavior and will foster further research since a desired 
simple tool is now available to design such nanoantennas 
with only little effort. It can be anticipated, for exam- 
ple, that the theory can be applied to design antennas 
supporting various resonances at predefined frequencies 
to respond on the desire to have multiresonant antennas 
at hand for applications in, e.g., Raman-sensors or ex- 
tremely broad-band resonators. The unique design tool 
of such circular patch nanoantennas consists in tailoring 
the dispersion relation and the complex reflection coef- 
ficient at will by carefully selecting a particular stack of 
layers. This provides a large degree of freedom which 
renders such optical nanoantennas unique. 
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APPENDIX 



Appendix A: Derivation of the Reflection Coefficient 



In this section, the complex reflection coefficient for Hankel-type surface plasmon polaritons (SPPs) for the given 
class of nanoantennas will be derived. For a definition of the geometry see again Fig.[2j At first, radially propagating 
plasmonic fields inside the resonator and a suitable ansatz for the outgoing field will be examined. At second, by 
calculating the Fourier components of the radiating field and by using the continuity of H " l , one can finally calculate 
the reflection coefficient integrating E™ ■ H™ at p = R. 



1. Plasmonic Modes and Outer Fields 

The studied nanoantennas obey a piecewise translational symmetry in z-direction. Thus, the structure consists of 
several layers which may be denoted by the subscript i in the following. It will be assumed that the antenna supports 
plasmonic surface modes propagating in radial direction. Hence, at least one of the materials is a metal. 

In each layer, one may split the fields into tangential (E^, Hi.t) and normal components (-Ej )Z , flj )Z ) relative to the 
interface of the discs. Then, the entire dynamic can be calculated using only the normal components. The tangential 
fields in the i-th layer are given by 

E,,t = jj [ki, z VtE i<z - unoHie z x VtH ix ] and 

k P 

H ,t = -rr [ki,zVtHi,z + ^e e l e 2 x V t Ei,z] (Al) 

kp 

with ki >z being the wave vector component in z-direction. In this formulation, all fields obey exp (ikt x z) and exp (— iu/t) 
dependencies^. 

The radial wave vector 



Kp — —^fJ-i^i K, 



%,z 



has to be independent of z since it is a conserved quantity in all layers. The various plasmonic modes of the antenna 
are thus characterized by different k p and it is convenient to define both 

fcspp = k p . 

In general, the structure may support several modes. Nevertheless, for the sake of simplicity, a further corresponding 
index will be suppressed. 

Furthermore, transverse-magnetic (TM) modes with H z = will be assumed and the materials shall be non- 
magnetic, i.e. jXi = 1. Thus, 

H z = , and 

H M = T2 — we o£ie z x VtEi,z 
fc SPP 



l 

2 — u)e eie z x 



fc S PP 



d p E iz e p + - {d v E uz )e 
p 



(A2) 



The only nonvanishing normal field component E z obeys the scalar two-dimensional Helmholtz equation 52 



A t E z + k 2 SPP E z = . 

Since fcspp is a constant throughout the layers, this differential equation generally applies to the corresponding 
plasmonic mode and the index i can be dropped. In rotational symmetry, the latter equation is known to have 
solutions in terms of Bessel and Neumann functions J m (kgppp) and N m (ksppp) modulus an angular variation of 
exp (irmp) which will be denoted by the index m in the following. 

For the calculation of the complex reflection coefficient one has to use Hankel functions which describe propagating 
rather than standing wave fields. They read as 

H U 2 (ksppp) = J m (ksppp) ± iN m (fcspp/o) 

and represent radially outgoing (1) and incoming (2) cylindrical waves in the given time-dependence. Since the fields 
are given in terms of E z it is convenient to introduce the reflection coefficient r m as 

E" 1 ^ (p, z) = A m (ksppp) ■ a (z) with 
A m (ksppp) = H m (k S p P p) + r rn ■ H m (k SP pp) 

where the boundary conditions at the antennas termination determine the "mode profile" a (z) in z-direction as well 
as the propagation constant /cspp- Furthermore, the superscript — was used to indicate fields for p < R. Likewise, + 
will be used in the outer region. 



From Eq. IA2I the angular component of the magnetic field is given by 



72^ + it 

= -, B m (k S ppp) e (z) a (z) 

kspp 

with 

B m (ksppp) = -, d p Am (ksppp) 

KSPP 

= DH m (k S ppp) + r m ■ DH 2 m (k SPP p) , 
DHimx)=d x HU 2 {x) . 

Having determined E™'~ and H" l: ~, a reasonable ansatz for the fields outside the antenna exhibiting the same 
symmetry properties, i.e. the same angular modulus m, can be constructed. For p > R, both fields under consideration 
are a superposition of outgoing Hankel-type SPPs whose amplitudes are given by certain Fourier coefficients c m (k z ). 
For each of these waves the wave vector in radial direction is given by 



K = \l edk l ~ k l ■ 

Here, a relative permittivity e d is assumed in the outer region. With this ansatz, the fields can be represented as 

ET' + (j>, z) = J" ^ (k z ) H m (Je d kl - B z p\ e ik * z dk z , 

/OO 
c m (k z ) ^-= (A4) 

-00 v erf — z 

■DH} n (\/e d kl - k*p) e lk > z dk z (A5) 

where the magnetic field follows from Eq. 



2. The Fourier Components of the Radiating Field 

In this section, the continuity of the magnetic field will be used to determine the Fourier components c m (k z ). Since 
B m , by definition, explicitly depend on r m , this will likewise hold for the c m (k z ). 
A Fourier transformation for p < R yields 

/OO 
H™>~ (p, z) e~ ik > z dz = -^B m (ksppp) ■ B- (k z ) 
-00 k S pp 

/OO 
e{z)a{z)e ±ik * z dz . 
CO 

The same operation for the outer region results in 

r H™>+ (p, z) e- ik **dz = i f c m (k) j^= DH m U e d k* - k*p) [°° e^^dzdk 



-00 



y/e d k% ~ k 2 

2irS{k-k z ) 



sje d kl - k 
Matching at p = R gives 



2^i • c m (k z ) j^$L= • DH m ( Je d kl kip 



-B m (k SP pp) B {k z ) = 27ri • c m (k z ) ■ 



^SPP V £ d fc - H 



■DHl Je d kl - k*R 
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further leading to the desired 

,. v _ & S ™ (k SPP R) ■ B- (k z ) 



B m (ksppR) _ y/e d k* - kjB- (k z ) 
27r-6 d -k S PP dh^ (^ Cdk 2 _ k 2jf 



3. Integration over E? ■ H™ 

The final step to calculate the reflection coefficient is the integration of the inner and outer E™ using H™'~ at 
p = R. Since there the factor B m (ksppR) is a constant, it is sufficient to use the field profile e (z) ■ a (z) alone. 
First, one finds 

/oo 
E" l ~ (R, z) e (z) a (z) dz = A m (k SPP R) ■ a with 
-OO 

a (z) • e (z) dz . 

> 

Performing the same integration over the mode field for p > R gives 

f E r z n '+ {R, z) e (z) a {z) dz = f c m (k z ) H m (je d k\ - k\R\ ■ f f e lk ' z t (z) a (z) dz) dk z 

= f Cm (k z ) H m Ue d k 2 - h? z R\ ■ B+ (k z ) dk z . 

and one can equate the last two results obtaining 

A m (k S ppR) ■ a = / c m (k z ) H m f y e d k% - k"\R J 
■B+{k z )dk z . 
Further using Eq. IA6l with the known c m (k z ), one finds 

A m (ksppR) ■ a = ■ I m (A7) 

2tt ■ e d ■ kspp 

with 



oo H m y/e d k% - klR 



— DH m y/e d k* - klR 



e d kl-klB-{k z )-B+{k z )dk z . 

Finally inserting the definitions of A m and B mi one can calculate the reflection coefficient as it is given in Eq. [5] in 
the main text, 

_ 2vr ■ e d ■ ksppa ■ H m {k S P P R) - DH m (k SPP R) ■ I m 
r,n -2tt • e d ■ ksppa ■ Hi (k SPP R) + DH m (k SPP R) ■ I m ' 

Also note that every rescaling in a (z) does not affect this result due to the linearity of the underlying theory. 

Appendix B: Comparison to Previous Results 

It is known that Hankel-type SPPs have the same dispersion relation as surface plasmons at planar surfaces^ which 
can be understood from the expansion of the Hankel functions for large arguments. 
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FIG. 10. The amplitude a) and phase b) of the reflection coefficient for different angular mode numbers m as a function of the 
radius R (solid lines). Assumed is a metal with a permittivity of e m = —15 + i0 + surrounded by air, i.e. td = 1. Results are 
compared to previous expressions that have been obtained for planar interfaces^ at normal incidence (dashed lines). 

The present theory can be applied to calculate the reflection coefficient of Hankel-type SPPs on top of a semi-infinite 
cylinder with the dielectric function e m (uj). The dispersion relation of a SPP at such an interface between a dielectric 
and a metal is then simply given byiiil 



fcspp(w) 




i(oj)ed 



(w) + e d 



(Bl) 



Since the expression is explicit, r m can also be explicitly derived, see App. [D] 

A further limiting case concerns the reflection of SPPs at an ordinary planar interface which is obtained for an 
infinite radius R — > oo. The result can be used to verify the theory in this limit since analytical expressions for such 
a specific geometry are available 3 ^. Results are outlined in Fig. [TUl There, both the modulus and the phase of the 
reflection coefficient are shown as a function of the radius of the semi- infinite cylinder for modes with different angular 
mode number m. It can be seen that the results as predicted by this theory converge for an infinite radius nicely 
towards the results previously reported 38 . 

Moreover, the theory discloses that there exists an intermediate region sustaining a very complex dependency with 
the emergence of a resonance behavior. It is associated with an abrupt decrease of the modulus and the emergence 
of a peak in the phase of the reflection coefficient. This region narrows for higher angular modes indicating a steep 
transition that might be of great practical use in the future, i.e. for a minor radius change the phase of the reflection 
coefficient changes abruptly; suggesting it for sensor applications. 

At a fixed radius, the reflection coefficient for discs of increasing thickness should converge to those of the semi- 
infinite disc. It is important to verify this point not only to check the consistency of the theory but also to see if one 
can approximately use the reflection coefficient of the semi-infinite case. 

Figure [TT1 shows that for d — > oo, the dipolar reflection coefficients of the even and odd modes approach that of the 
semi-infinite approximation for a radius of k$R = 20 and e m = — 15 + i0 + , e^ = 1. 

The absolute value of r m converges more rapidly in the given case. However, one can see that already at ko ■ d > 3, 
the values do not vary strongly with respect to the semi-infinite value which would justify to use this approximation 
here. 

From all these considerations it may be concluded that in limiting cases the theory reproduces previous results; 
being an indication for its performance. It can be also seen that although the theory might have restrictiions in its 
formulation; it turns out to be rather general since many other geometries can be deduced from it while considering 
the limits of the theory. 



Appendix C: Explicit Plasmonic Modes 

It was shown 39 that a Hankel-type plasmonic mode at a dielectric-metal surface has the same dispersion relation 
as the one known for a plane-wave plasmonic mode. This is the general case as can be seen by the infinite-radius 
approximations of the cylindrical Hankel functions. In this section, the analytic expressions for plasmonic modes as 
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FIG. 11. Comparison of ri for e m = — 15 + i0 + and e<j = 1 with the semi-infinite cylinder model at koR = 20 which is approached 
as the thickness d goes to infinity both in absolute value and phase, see a) and b). 

used in the main body of the manuscript will be explicitly given for convenience - for semi-infinite cylinders, discs 
with finite thicknesses and a stacked antenna consisting of two metallic discs separated by a dielectric spacer. 



1. Semi-Infinite Cylinder 



Assumed is a cylinder of radius R terminated at z = 0. Its relative permittivity is then described by 

y e m z < 
where 

e' d > and 
e m < ~ e d ■ 
In this case, a Hankel-type surface plasmon polariton exists with the mode profile 

' ' e m e~ ldZ z>0 



a (z) = 



e d e^ z z < ' 



e (z) a (z) = e d e r , 



e ~i d z z > o 
e^ z z<0 



and the well-known dispersion relatio 



6.7 



Kgpp 



— — k n and 



td 



■y 2 - k 2 -fk 2 
-k 2 



(d + ^r 



(CI) 



By definition the 7d/ m are related to penetration depths of the plasmonic fields d pm ^ — ( 7 m / d 



2. Disc with Finite Thickness 



In this subsection, the case of a metallic disc with radius R surrounded by another dielectric material will be 
considered. Then, for p < R, the relative permittivity is given by 



e(z) 



e d z> d/2 

e m \z\ < d/2 
e d z < -d/2 
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In this case, there are two transverse magnetic modes supported by the system. They can be described by 

e m exp [-7d(z- f)] z > d 
a± (z) = I Te™ exp [j d (z + f )] z < 

and the dispersion relation is given by 



hyp T ( 7mZ ) , , rf 



tanh ( 7 m - 



±1 
7™ £d 



7^ 



with 



2 _ i 2 .2 

7t — %PP — e i K 



Modes which are odd in E z , thus odd in a(z) are called even and denoted here by a + and &+. In this case, 
hyp_ (x) = sinh(x). An illustration of this mode can be found in Fig. [TJ]a). For odd modes, the corresponding 
conventions hold along with hyp + (x) = cosh(x). Then, one finds a and B~ (k) ■ B + (fc) by direct computation. 

3. Stacks with Two Metallic Discs 

Assumed is that two metallic discs with radius R and thickness d m are separated by some dielectric with thickness 
d d . Furthermore, this simple stack is surrounded by the same dielectric, hence 

J e d z < 0, d m < z < d m + d d and 2d m + d d < z , 
e (z) = < 

[e m < z < d m and d m + d d < z < 2d m + d d 

The structure supports several Bessel-type plasmonic modes that can be characterized by 

'e m Ci +ti e^ z z<0 

e d C 2 +, l e^^ z - d ^ + e d C 2 -,ie-^ iZ < z < d m 

at (z) = < e m C 3+ ^^ z ~ d ^ d ^ + e m C 3 -, ie -^ z - d ^ d m <z<d m + d d 

e d C 4+ ^^ z - 2d ™- d ^ + e d C 4 -, l e-^* ( - z - d ™- d ^ d m + d d < z < 2d m + d d 

l e m C 5 _, 4 e-'w-'(*- M »- , ''0 2d m + d d < z 

with r i^ n i di = k 2 SPPi — e m /d,iko- The dispersion relation and the coefficients Ci±^ can be found by applying the 

transfer matrix method 6 which can also be used for more complex stack compositions. 

In the study of Section IIVC1 a symmetric situation was used with d d — d m such that only two dominating modes 
have to be considered. The first one is characterized by 

Ci+,i = Cs-,1' 03+,! = U3- t i, C2-,i = C*4+,i and C'2+,1 = C±- t i . 

With the explanations given in the preceding subsection, it is natural to call this mode antisymmetric since it is 
symmetric in E z ex a(z). The second mode is given due to the relations 

Cl + ,2 = — C^-,2, C3+,2 = —G'3-,2, G'2-,2 = —Ci+.2 and C2+,2 = — Ck-,2 

and consequentially termed symmetric. The modes are illustrated in Fig. [T^b); see again Fig. [5] for an illustration of 
the nanoantenna. However, the explicit coefficients C;±^ and corresponding dispersion relations for kspp^i turn out 
to be rather complicated to be presented here. 

Appendix D: Explicit Reflection Coefficient for a Semi-Infinite Cylinder 

The underlying physical situation of SPPs on a semi-infinite cylinder is a limit case in two aspects: For large radii, 
the reflection coefficient is converging towards that of SPPs for normal incidence at a planar interface^. Furthermore, 
it is equivalent to the reflection of Hankel-type SPPs on cylinders of large thicknesses. This relation was used in the 
main body of the manuscript to verify the calculations. 

In the given case, the dispersion relation takes an explicit form, see Eq. IC11 and it is possible to find a more explicit 
formula for the reflection coefficient which will be derived in the following. 
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FIG. 12. Modes that could be observed in the numerical experiments performed. The material of the metals is taken to be 
silver with e m w —8.8 + 0.3i at v — 625 THz. In a), one can see e(z)a(z) for the even mode of a single disc with 30 nm 
thickness. The same quantity is outlined in b) for the symmetric and antisymmetric modes of a nanoantenna consisting of two 
silver discs with 30 nm thickness and also 30 nm air cavity in between. These modes can be understood as symmetric and 
antisymmetric coupling between single discs obeying an even mode. One can directly see that in this case the antisymmetric 
mode is strongly confined inside the cavity. On the contrary, E z of the symmetric mode is mostly located outside the metal 
discs. 



As a preliminary step, one has to calculate B (fc) • B + (fc) and a using the permittivity e (z) and mode profile a (z) 



B- (k) ■ B + (k) = e 2 d e 



2,2 



{id + 1m) 



m {k 2 +l 2 d )(k*+il) ' 



Zldln 

2 2 
e d e m 



(idtd + 7m£m) 

1 1 

+ 



2 [ Cmlm tdld 



Taking Eq. IA7l from the derivation of the complex reflection coefficient, one can further simplify 



A m {ksppR) 



^d^m 



^mlri 



Cdld 



Bm (ksppR) 
2-k ■ e d ■ kspp 



e dk k z € d € m 



m 



e d k 2 - k\R 



DH m (y/e d k% - klR 

\2 



(ld+lrn) 



{k 2 +l 2 d )(k 2 +i 2 m ) 



dk. 



to 



An (ksppR) = 



B m (k SP pR) {ld + lm) 



■K 


£d ■ kspp 


c, 


1 




\/ e d k - 


JU2 





1 

tdld 



(k 2 + l 2 d )(k 2 +i 2 m ) 



dk z 



H m (y/e d k 2 ~k 2 R 



DH m ^e d k 2 - klR) 



Using 



{lm + Id) 

tdid 



e™7 



V e d^m ' kgpp ■ 



ed 



£d 



which follows directly from the definitions of the ji, one obtains 



a n m ® m ( k SPpR) i 

Am {ksppK) = ■ V- £ d e r 



,2 ^ m ^d 

■ Kgpp ■ 



H m (Vtdk'Z - k*R 



e m + £d J- 



oc DH l 



V e dk k z 

(P+7 2 )(fc 2 +7 2 n) ^ 

_ B m (k S ppB) 



Q7T 



' %/-£d£r 



^7n^d 



,00 H, 



e m + e d e m + e d 



'e d kl - klR) 



I ( y/e d k% - klR 



*> DHl, ( v^eF^? 



and further 



using the abbreviations 






7rV-e d e m An (k SPP R) = S m (k SPP R) ■ e 2 p ■ ( 1 — j 1^ 



/"DC 



rOO 



v^ 



gp 1 (VJI^H^fl) | 



e„ = 



** —00 






6d + £r 



Then, the reflection coefficient is given by 

Ti^e d e 



1 ID 



• H} n {k SPP R) - DHl (ksppR) e 2 p -(l- ^) 7= 

• HI (ksppR) + DHl (ksppR) -4-(l-^y 
ies to the reflection of SPPs at a planar interface 38 . 



-*V=^ ■ ff£ (k SPP R) + DHl (ksppR) ■ 4 ■ (l - ij) /£ 
This result has, as expected, similariti ' " '" ' '•"*" ; ' ' 
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